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ABSTRACT 

The  interference  reduction  capability  of  antenna  arrays  and  the  power  control  algorithms 
have  been  considered  separately  as  means  to  increase  the  capacity  in  wireless  communi¬ 
cation  networks.  The  MVDR  (Minimum  Variance  Distortionless  Response)  beamformer 
maximizes  the  Carrier  to  Interference  Ratio  (CIR)  when  it  is  employed  in  the  receiver  of  a 
wireless  link.  In  a  system  with  omnidirectional  antennas,  power  control  algorithms  are  used 
to  maximize  CIR  as  well.  In  this  paper  we  consider  a  system  with  beamforming  capabilities 
in  the  receiver,  and  power  control.  An  iterative  algorithm  is  proposed  to  jointly  update 
the  transmission  powers  and  the  beamformer  weights  so  that  thet  converge  to  the  jointly 
optimal  beamforming  and  transmission  power  vector.  The  algorithm  is  distributed  and  uses 
only  local  interference  measurements.  In  an  uplink  transmission  scenario  it  is  shown  how 
base  assignment  can  be  incorporated  in  addition  to  beamforming  and  power  control  such 
that  a  globaly  optimum  solution  is  obtained.  The  network  capacity  increase  and  the  saving 
in  mobile  power  achieved  by  beamforming  are  evaluated  through  numerical  study. 

Keywords:  Adaptive  Beamforming,  Power  Control,  Space  Division  Multiple  Access  (SDMA). 


This  work  is  supported  in  part  by  the  NSF  grant  MIP9457397,  NSF  grant  NCR-9406415,  CAREER  award 
NCR-9502614,  and  AFOSR  95-1-0061. 


I.  Introduction 


Cochannel  interference  is  one  of  the  main  impairments  that  degrades  the  performance  of  a 
wireless  link.  Power  control  and  antenna  beamforming  are  two  approaches  for  improving  the 
performance  in  wireless  networks  by  appropriately  controlling  the  cochannel  interference. 

In  power  control,  the  transmitter  powers  are  constantly  adjusted  by  increasing  if  the  Signal 
to  Noise  Ratio  (SNR)  is  low  and  decreasing  if  the  SNR  is  high,  so  that  the  quality  of  weak 
links  is  improved.  Receivers  employing  antenna  arrays  may  adjust  their  beam  patterns  such 
that  they  have  maximum  gain  towards  the  directions  of  their  transmitters  and  minimum 
gain  towards  the  other  directions  so  that  the  aggregate  interference  power  is  minimized. 

Previous  work,  discussed  later  in  more  detail,  addressed  separately  the  problems  of  power 
control  for  optimal  interference  balancing  and  beamforming  to  maximize  the  SNR.  In  this 
paper  we  consider  the  joint  problem.  We  consider  a  set  of  cochannel  links,  the  receivers 
of  which  employ  antenna  arrays,  e.g.,  a  set  of  cochannel  uplinks  in  a  cellular  network.  An 
algorithm  is  provided  for  computing  the  transmission  powers  and  the  beamforming  weight 
vectors,  such  that  a  target  SNR  is  achieved  for  each  link  (if  it  is  achievable)  with  minimal 
total  transmission  power.  The  algorithm  is  decentralized  and  amenable  to  a  distributed 
implementation.  It  operates  as  follows.  For  a  fixed  power  allocation,  each  base  station 
maximizes  the  CIR  using  MVDR  beamformer.  For  this  fixed  beamforming  pattern,  the 
mobile  powers  are  updated  to  reduce  the  cochannel  interference.  This  operation  is  done 
iteratively  until  both  the  vector  of  transmitter  powers  and  the  weight  coefficients  of  the 
beamformers  converge  to  the  jointly  optimal  value.  In  the  case  that  each  transmitter  has 
multiple  options  regarding  which  receiver  to  engage,  e.g.  each  mobile  may  select  its  base 
station  among  a  set  of  possible  options,  the  algorithm  easily  extends  to  find  the  joint  optimum 
power,  base  station  and  beamforming. 

The  application  of  antenna  arrays  has  been  proposed  in  [1]  to  increase  network  capacity  in 
CDMA  systems.  This  work  was  based  on  equal  received  power  from  all  users  in  base  stations. 
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In  [2]- [5]  centralized  power  control  schemes  have  been  proposed  to  balance  the  Carrier  to 
Interference  Ratio  (CIR)  or  maximize  the  minimum  CIR  in  all  links.  These  algorithms  need 
global  information  about  all  link  gains  and  powers.  The  distributed  power  control  algorithm 
which  uses  only  local  measurements  of  CIR  was  presented  in  [6]- [9].  In  [10]  and  [11]  the 
combined  base  station  assignment  and  power  allocation  was  used  to  increase  uplink  capacity 
in  wireless  communication  networks.  In  those  works  it  was  shown  that  if  there  exists  at  least 
one  feasible  base  station  assignment,  the  proposed  algorithms  will  converge  to  the  optimal 
base  station  assignment  and  power  allocation  in  the  sense  that  it  minimizes  the  interference. 
Zander  has  noted  [4]  in  some  situations  the  minimum  protection  ratio  is  not  achievable  by 
all  links.  In  other  words,  there  is  no  feasible  solution  to  the  power  control  problem.  In  [4] 
he  suggested  removing  a  minimum  number  of  links  to  enable  the  other  links  to  establish 
acceptable  connections. 

The  paper  is  organized  as  follows.  In  Section  ,  we  describe  the  system  model  and  existing 
power  control  algorithms.  Section  considers  the  beamforming  problem  in  a  network  of  users. 
In  Section  ,  we  consider  power  allocation  and  beamforming  as  a  joint  problem.  An  iterative 
algorithm  for  the  joint  problem  is  presented  which  converges  to  the  optimal  solution,  i.e., 
the  allocated  powers  are  minimum  among  all  sets  of  feasible  power  allocations.  In  Section 
,  we  consider  an  uplink  scenario  where  in  addition  to  power  control  and  beamforming  base 
station  assignment  is  performed  as  well. 

In  Section  a  simulation  study  is  done.  We  will  show  that  our  method  increases  the 
capacity  in  cases  when  users  are  uniformly  dispersed  around  the  network,  and  when  the 
users  are  added  in  a  locally  congested  area  in  the  network.  The  simulation  results  show  that 
using  antenna  arrays  with  four  elements,  the  speed  of  convergence  is  increased  up  to  three 
times  and  the  capacity  is  increased  up  to  five  times  compared  to  the  case  of  joint  power 
control  and  base  station  assignment  with  omnidirectional  antennas. 
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II.  System  Model  and  Power  Control  Problem 


A  set  of  M  transmitter- receiver  pairs  which  share  the  same  channel  is  considered.  The 
shared  channel  could  be  a  frequency  band  in  FDMA,  a  time  slot  in  TDMA,  or  even  CDMA 
spreading  codes.  The  link  gain  between  transmitter  i  and  receiver  j  is  denoted  by  Gij,  and 
the  zth  transmitter  power  by  Pi.  For  an  isotropic  antenna  with  unity  gain  in  all  directions, 
the  signal  power  received  at  receiver  j  from  transmitter  i  is  GijPi,  as  illustrated  in  Fig.  1. 
It  is  assumed  that  transmitter  i  communicates  with  receiver  i.  Hence  the  desired  signal  at 
receiver  i  is  equal  to  GuPi ,  while  the  interfering  signal  power  from  other  transmitters  to 
receiver  i  is  I  =  Yjjfr  GjiPj-  The  CIR  at  the  ith  receiver  is  given  by 


r,= 


GuPi 


GjiPj 


If  we  neglect  thermal  noise,  the  quality  of  the  transmission  link  from  transmitter  i  to  receiver 
j  depends  solely  on  Tj.  The  quality  is  acceptable  if  Tj  is  above  a  certain  threshold  70,  the 
system  protection  ratio,  which  is  determined  based  on  the  signaling  scheme  and  the  link 
quality  requirements  (target  bit  error  rate).  Hence  for  acceptable  link  quality, 


G  P 

V-J n 1  1 


GjiPj 

In  matrix  form,  (1)  can  be  written  as  follows, 


>  7o- 


(1) 


1_ 

7o 


P  >  F  P, 


(2) 


where  P  =  [P\,  P2, . . . ,  Pm}T,  is  the  power  vector,  and  F  is  a  non-negative  matrix  defined  as 


Fij  =  < 


0  if  j  =  i 

>  0  if  j±i 

The  objective  of  a  power  control  scheme  is  to  maintain  the  link  quality  by  keeping  the  CIR 
above  the  threshold  70,  that  is,  to  adjust  the  power  vector  P  such  that  (2)  is  satisfied.  This 
problem  has  been  studied  extensively  recently  [2]-[ll].  Given  that  F  is  irreducible,  it  is 
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known  by  Perron-Frobenius  theorem,  that  the  maximum  value  of  y0  for  which  there  exists 
a  positive  P  such  that  (2)  is  satisfied  is  ^y,  where  p(F)  is  the  spectral  radius  of  F  [12]. 
According  to  this  theorem,  the  power  vector  that  satisfies  (2)  is  the  eigenvector  corresponding 
to  p(F)  and  is  positive.  Considering  the  thermal  noise  at  the  receivers,  the  CIR  at  the  Ah 
receiver  is  modified  as 

P  _ _ GuPi _ 

4  ~  E&GjiPj  +  IV,  ’ 

where  jV*  is  noise  power  at  the  Ah  receiver.  The  requirement  for  acceptable  link  quality  is 
again 

R  >  7o,  1  <  i  <  M, 


or  in  matrix  form 


[I  -  7qF]P  >  u, 


(3) 


where  I  is  an  M  x  M  identity  matrix,  and  u  is  an  element-wise  positive  vector  with  elements 
Ui  defined  as 


Ui  — 


7o  Ni 


Gu  ’  - 


1  <i<  M. 


The  CIR  threshold  70  is  achievable  if  there  exists  at  least  one  solution  vector  P  that  satisfies 
(3).  The  power  control  problem  now  is  as  follows: 


minimize  ^  P,, 


subject  to  [I  —  7oP]P  >  u. 

It  can  be  shown  that  if  the  spectral  radius  of  F  is  less  than  A,  the  matrix  I  —  70 F  is  invertible 
and  positive  [12].  In  this  case  the  power  vector 

P=:[/-7oP]-1u  (4) 


solves  the  optimization  problem. 
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A  centralized  power  control  algorithm  ([4] -[5])  solves  (4)  by  requiring  all  link  gains  in  the 
network  and  noise  levels  at  receivers.  In  [6]- [8]  a  decentralized  solution  to  the  power  control 
problem  is  proposed  that  solves  (4)  by  performing  the  iterations 

Pn+1  =  7oFPn  +  u.  (5) 

The  right  hand  side  of  (5)  is  a  function  of  the  interference  at  the  ith  receiver,  and  the  link 
gain  between  each  receiver  and  its  transmitter.  That  is,  there  is  no  need  to  know  all  the 
existing  path  gains  and  transmitter  powers  in  order  to  update  the  powers.  At  each  iteration, 
transmitters  update  their  powers  based  on  the  interference  measured  at  the  receivers  and 
the  link  gain  between  each  transmitter  and  its  own  receiver,  which  can  be  measured  from 
the  information  sent  in  the  control  channel.  It  has  been  shown  in  [6]-[8]  that  starting  from 
any  arbitrary  power  vector,  this  solution  converges  to  the  optimal  solution  P. 

III.  Antenna  Array  and  Beamforming 

An  antenna  array  consists  of  a  set  of  antennas,  designed  to  receive  signals  radiating  from 
some  specific  directions  and  attenuate  signals  radiating  from  other  directions  of  no  interest. 
The  outputs  of  array  elements  are  weighted  and  added  by  a  beainformer  as  shown  in  Fig.  2 
to  produce  a  directed  main  beam  and  adjustable  nulls.  In  order  to  reject  the  interferences, 
the  beamformer  has  to  place  its  nulls  in  the  directions  of  sources  of  interference,  and  steer  to 
the  direction  of  the  target  signal  by  maintaining  constant  gain  at  this  direction.  A  sample 
antenna  array  pattern,  which  is  depicted  in  Fig.  3  shows  this  effect. 

Assume  the  array  consists  of  K  +  1  elements,  and  its  response  at  the  direction  of  arrival  9 
is  a (9)  =  [ao(0),  ai(6),  ...  ,  clk(9)],  where  a,k{9)  is  the  response  of  the  A;th  antenna  element  to 
the  signal  s(t)  from  direction  9.  The  signal  vector  received  due  to  s(t)  can  be  written  as 

x(f)  =  a  (9)VPs(t), 
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where  P  is  the  received  power  of  the  source.  The  (K  + 1)  x  1  vector  a(0)  is  called  the  spatial 
signature  or  array  response  for  the  source  at  direction  9. 

Now  consider  a  cochannel  set  consisting  of  M  transmitter  and  receiver  pairs.  Denote  the 
ith  array  response  to  the  jth  source  by  an.  The  received  signal  at  the  ith  receiver  can  be 
written  as  follows 

M 

=  E  V^siWati  + 

3= 1 

where  n*(f)  is  the  thermal  noise  vector  at  the  input  of  antenna  array  at  the  ith  receiver,  and 
Pj  is  the  power  of  the  jth  transmitter.  Consider  the  problem  of  beamforming  as  to  maximize 
CIR  for  a  specific  link,  which  is  equivalent  to  minimizing  the  interference  at  the  receiver  of 
that  link.  In  order  to  minimize  the  interference,  we  minimize  the  variance  or  average  power 
at  the  output  of  the  bearnformer  subject  to  maintaining  unity  gain  at  the  direction  of  the 
desired  signal.  We  can  write  the  output  of  the  bearnformer  at  the  ith  receiver  as 

e*  =  wf  Xj, 

where  w*  and  x*  are  the  beamforming  weight  vector  and  the  received  signal  vector  at  the 
ith  receiver,  respectively.  The  average  output  power  is  given  by 

Si  =  E{wfxjxfwj}  =  wf  E{xjxf  }wi  =  wf$jWj,  (6) 

where  <f>j  is  the  correlation  matrix  of  the  received  vector  Xj.  If  the  message  signals  Sj(t)  are 
orthogonal  to  each  other,  the  correlation  matrix  is  given  by  [13] 

=  E  pjGjiajiafi  +  Ni 1  +  PiGu&ii&fi 

=  in  T  RjGjjajjajj ,  (7) 

where 

=  E  PjG3ia3iafi  +  NG 
j¥=i 
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is  the  correlation  matrix  of  unwanted  signals,  and  Nt  is  the  noise  power  at  the  input  of  each 
array  element.  Combining  (6)  and  (7),  we  obtain  the  received  signal  plus  interference  power 
as  a  function  of  weight  vector  wp 

Si  =  PiGu  +  Niwfwi  +  Y  PjGjiw?a.jiafiWi.  (8) 

Here  we  use  the  fact  that  the  gain  at  the  direction  of  interest  is  unity,  i.e,  =  1.  The 

first  term  in  (8)  is  the  received  power  from  the  signal  of  interest,  while  the  other  terms  are 
related  to  the  interference  and  noise.  Therefore  the  total  interference  is  given  by 

h  =  YGjiwfajiafiwiPj  +  Ni  wfwi- 

The  goal  of  beamforming  is  to  find  a  weight  vector  w*  that  minimizes  the  interference 
subject  to  vrfaa  =  1.  It  can  be  shown  that  the  solution  is  [14] 

<I>— ^a-  ■ 

-  _  ^ i  aii 

Wi  ~  affsrw 

The  antenna  gain  at  the  direction  of  interest  is  unity,  that  is,  the  received  desired  signal  is 
unaffected  by  beamforming.  Therefore,  the  CIR  at  the  ith  receiver  is  given  by 

Tj  =  PlGtlafl^all. 

In  order  to  calculate  the  received  power  for  transmitter  j,  we  have  to  multiply  the  trans¬ 
mitter  power  by  the  antenna  power  gain  in  addition  to  the  propagation  path  gain,  i.e., 


G jiG  a%  (wj,  aji), 


where  Ga,  (w»,  a^)  =  wf  a^a^w p  Then  the  maximum  CIR  at  the  zth  receiver  can  be  written 
as 


Ti 


PiGu 


GjiGai  (wp  aji)Pj  +  Ni  wfw' 


(10) 
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Here  we  assumed  that  the  array  response  to  the  source  of  interest  is  known.  The  array 
response  is  obtained  by  the  estimation  of  the  direction  of  arrival  (DOA).  But  in  wireless  net¬ 
works  usually  the  number  of  cochannels  is  much  larger  than  the  number  of  array  elements. 
Therefore,  conventional  DOA  estimation  methods  like  ESPRIT  and  MUSIC  are  not  appli¬ 
cable.  However,  there  exist  several  schemes  that  can  be  used  to  estimate  the  array  response 
in  TDMA  ([15]),  FDMA  ([16]),  and  CDMA  networks  ([13]),  without  estimating  the  DOA. 


IV.  Jointly  Optimal  Power  Control  and  Beamforming 


In  the  joint  power  control  and  beamforming  problem  the  objective  is  to  find  the  optimal 
weight  vector  and  power  allocations  such  that  the  CIR  threshold  is  achieved  by  all  links, 
while  each  transmitter  keeps  the  transmission  power  at  the  minimum  required  level  to  reduce 
the  interference  to  other  users.  The  CIR  at  the  ith  receiver  is  given  by 


r,;  = 


PiGii 


Ejfr  GjiGa,  (w i,  a ji)Pj  +  Ni  wf  Wj 
The  optimization  problem  now  is  defined  as 


mm 

W.P 


M 

E* 

i= 1 


subject  to  R  >  7o,  (i  =  1, 2, . . . ,  M). 


(11) 


This  constraint  can  be  presented  in  matrix  form  as: 


[/  -  l0Fw}P  >  u, 


where 


[FWh  = 


0  if  j  =  i 

G^Gafv/^a^)  >  q  0^erwjse 


and  u  is  an  element-wise  positive  vector  with  elements  Ui  defined  as 

Awfw, 


Ui  = 


'l  ,yi 

Gu 


(*  =  1,2, ...,M). 


8 


Assume  that  there  is  a  set  of  weight  vectors  W  =  {wx,  W2, . . . ,  w m},  for  which  p(Fw )  <  A. 
The  matrix  I  —  70 Fw  is  then  invertible  and  =  [I  —  70 Fw]~1  u  minimizes  the  objective 
function  in  the  optimization  problem  for  the  fixed  weight  vector  set  W.  For  any  feasible 
pair  (W,P„,),  the  vector  Pw,  can  be  computed  as  the  limit  of  the  following  iteration 


nn+1  _  .  V-  GkiGat(wi,aki)  nn  7oATiwiHwi 
F  =7o  2^ - FT- - F  + 


k^i 


Gu 


Gu 


(i  =  l,2,...,M). 


(12) 


The  above  iteration  is  the  same  as  the  distributed  power  control  algorithm  (see  [6]- [9]),  in 
which  the  path  gain  Gki  is  replaced  by  the  multiplication  of  the  path  gain  and  antenna  gain, 
and  the  noise  power  is  replaced  by  the  weighted  sum  of  the  noise  powers  at  the  input  of  each 
antenna  array.  Denote  the  iteration  in  (12)  as 


P”  =mw(Pn_1). 


Starting  from  any  initial  power  vector  P°,  the  mapping  mw  will  converge  to  the  optimal  power 
vector  Pw  which  is  the  fixed  point  of  the  mapping,  i.e.,  limn_>oo  Pn  =  P^,  P^  =  mw(Pw).  The 
objective  in  joint  optimal  beamforming  and  power  control  problem  is  to  find  the  beamforming 
set  W  among  all  feasible  beamforming  sets,  in  such  a  way  that  Pw  is  minimal. 

The  level  of  cochannel  interference  depends  both  on  the  gain  between  interfering  transmit¬ 
ters  and  receivers,  as  well  as  on  the  level  of  transmitter  powers,  i.e.,  the  optimal  beamforming 
vector  may  vary  for  different  powers.  Hence,  beamforming  and  power  control  should  be  con¬ 
sidered  jointly. 

Let  (W,  P)  be  the  weight  vector  and  power  vector  pair  which  achieves  the  minimum  in 
(11).  In  the  sequel,  we  present  an  iterative  algorithm  for  adjusting  P  and  W  simultaneously, 
and  we  will  show  that,  starting  from  any  arbitrary  power  vector,  it  converges  to  the  optimal 
solution  (W,  P).  The  iteration  step  for  obtaining  (Pn+1,  Wn+1)  given  (Pn,  W”)  is  as  follows: 
Algorithm  A: 

1.  wf+1  is  computed  at  each  receiver  i  such  that  the  cochannel  interference  is  minimized 
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under  the  condition  of  constant  gain  for  the  direction  of  interest,  i.e., 


w”+1  =  arg  min{£  G3lGa%  (w* ,  a,*) P?  +  N{  wfw * } ,  (z  —  1,  2, ,  M), 

subject  to  WjHajj  =  1, 

where  Pn  is  the  power  vector  updated  at  step  n. 

2.  The  updated  power  vector,  Pn+1,  is  then  obtained  by 

pn+i _ _  v  ^(w^.ajti)  7o-W(w”+1)ffw"+1 


k^i 


Ga 


._7X  — J—  1 


that  is,  by  performing  one  iteration  with  the  mapping  mw*  on  the  power  vector  Pn. 
Combining  two  iteration  steps  in  the  algorithm,  we  obtain  the  power  vector  update  in  a 
single  step: 


Prl  =  min{T»EG“G"iW,'a“)f?  + 


W*  ‘"i &  G> 


G, 


subject  to  Wtf/a a  =  1. 


(13) 


We  denote  the  iteration  in  (13)  as 

Pn+1  =  m(P"). 

Theorem  1:  The  sequence  (Pn,Wn),  (n  =  1,2,.. .)  produced  by  the  iteration  (13),  starting 
from  an  arbitrary  pair  (P°,  W°),  converges  to  the  optimal  pair  (P,  W). 

In  order  to  prove  theorem  1,  it  is  shown  first  that  the  theorem  holds  when  the  iteration  starts 
from  the  power  vector  P°  =  0. 

Lemma  1:  For  any  two  power  vectors  P^  and  P#  such  that  P^  <  P#  the  following  holds: 

(a)  m(PA)  <  mw(PA),  Vw 

(b)  mw(PA)  <  mw( PB),  Vw 

(c)  m(PA)  <  m{P B) 
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Proof:  (b)  and  (c)  can  be  concluded  immediately  from  the  fact  that  the  coefficients  in  the 
mappings  m  and  mw  are  positive,  (a)  holds  since  in  the  mapping  m,  we  are  minimizing  the 
power  vector  P  over  all  possible  weight  vectors  w. 

□ 

Theorem  2:  The  sequence  Pn  by  iteration  (13)  and  initial  condition  P°  =  0,  converges  to 
the  fixed  point  of  the  mapping  m ,  P. 

Proof:  We  define  two  power  vector  sequences  Pn  and  PJ*  produced  by  the  mappings  m  and 
mw,  respectively,  with  zero  initial  condition.  That  is, 


P”+1  =  m(Pn),  P°  =  0, 


and 


PS+1  =  m*(PS), 


=  0. 


The  power  vector  sequence  Pn  is  non-decreasing.  In  order  to  show  this  we  observe  that 
P1  =  m(P°)  =  u  >  0,  i.e.,  P°  <  P1,  and  if  Pn_1  <  Pn  lemma  1(c)  implies  m(Pn_1)  <  m(Pn) 
or  P"  <  Pn+1.  By  induction  we  conclude  that  Pn  is  a  non-decreasing  sequence.  We  can 
follow  the  same  steps  to  prove  that  the  sequence  P£  is  also  non-decreasing. 

We  start  the  mappings  m  and  mw  from  the  same  starting  vector  P°  =  P£,  =  0.  By  lemma 
1(a)  m(P°)  <  m*( P^)  or  P1  <  Pi,  and  if  Pn  <  P£,  by  lemma  1(a), (b)  m( P")  <  mw{J?r~!) 
or  Pn+1  <  P(f+1  for  all  n.  That  is,  by  induction  we  may  write  P"  <  P^,,  (for  n  =  1,2,...). 
Since  W  is  the  optimal  beamforming  set,  the  sequence  P£  will  converge  to  the  optimal  power 
vector  P,  i.e., 


lim  PS,  =  P. 


Hence  Pn  is  a  non-decreasing  sequence  and  it  is  bounded  from  above  by  P,  so  it  has  a  limit 
denoted  by  P*.  It  is  shown  in  the  following  that  P*  =  P.  Since  Pn  approaches  P*  from 
below  we  can  write  Pn  =  P*  +  en,  where  en  =  P"  -  P*  is  non-positive  and  approaches  zero. 
Now  we  show  P*  is  the  fixed  point  of  the  mapping  m.  Since  e  <  0,  Pn+1  =  m(Pn)  implies: 
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pr  =  ^EG-^S^W  +  4)  +  2 »1} 


w  •  * — *  (7.. 


.  ^GwGai(wi,a«)  ToMwfwi  ^  GkiGai(-whaki) 
=  rmn{7o  E - 75 - Pk  + - 7^ - I-  7o  E  ~ 


Wi  G- 


k^=i  Gii 


4) 


.  .  (  v- GfciGai(wi,afci)  To^wfwj 

<  mm  {70  - — - Pfe  + - ^ - } 


•Wi  '  ■  -  t—*  Q.. 

k^i  ^ll 


Ga 


(14) 


Also  Pn+1  =  P*  +  en+1,  concludes  that  P*  +  en+1  <  m(P*). 
Since  en+1  approaches  zero, 


P*  <  m(P*). 


(15) 


Again  from  Pn+1  =  m(Pn),  we  conclude: 


PI 


n+ 1  _ 


min 

w 


rin(7o  E 


GkiGa.  (wii  &ki)  /  D* 


k^i 


G, 


(pk  + «?)  + 


7oAiwfw; 

G« 


} 


^  .  r  ^  GwG0i(wi,aJM)  m  ,  7oA/'»wfwi  r  ^  GfciGfl,(wi,aw)  n1 

>  min {7o  E - - P*  +  “ — 75^" —  I  +  min{7o  E - — ~ek > 


k^i 

>  rrii(P*)  +  nii(en), 


G, 


k^ti 


(16) 


where 


Consequently, 


rhi(en) 


min{7oE 


GkiGai  (wj,  a**) 


G« 


P*  +  en+1  >m(p*)  +  m{en). 
m(en)  and  en  both  approach  zero,  therefore 


P*  >  m(P*). 


(17) 


12 


Finally  from  (15)  and  (17),  P*  =  m(P*). 
Let 


wf  =  arg  min{^  G3lGai  (w* ,  a^) P*  +  Nt  wf w* } ,  (*  =  1 , 2, . . . ,  M) , 


subject  to  =  1, 


and 


D*  r  V-  GkiGati-WuHbi)  D*  ,  70-Wi(Wi)HWi 

W  -  \7o  2_,  - 7i - P*  + . 7? - }• 

A:#  U*» 

By  definition  w*  =  Wj  and  consequently  P*  =  P.  That  is,  the  sequence  P"  converges 
to  the  optimal  power  vector  P.  Furthermore,  the  weight  vectors  converge  to  the  optimal 
beamforming. 

□ 

Proof  of  theorem  1:  Now  we  will  show  that  a  power  vector  sequence  starting  from  any 
initial  power  vector  converges  to  the  optimal  power  vector  P.  We  consider  the  sequence 
pn+i  _  with  the  arbitrary  initial  power  vector  P°. 

Assume  there  exists  a  feasible  pair  (P,W).  The  power  vector  iteration  for  this  pair  is 
given  by 

PJ+1  =  m”(P”),  (n  =  0,1,...).  (18) 

The  optimality  of  W  implies  that  lim^oo  P£  =  P,  where  P  is  the  fixed  point  of  the  mapping 
defined  in  (18).  Assume  that  both  sequences  start  from  the  same  point,  i.e.,  P°  =  P^. 
Lemma  1(c)  implies  m(P°)  <  mw(P%)  or  P1  <  Pi.  If  Pn  <  P£,  then  m(Pn)  <  mw( P(~)  or 
pn+i  <  p^+1.  Hence  by  induction  we  have 

P”<P£,  (n  =  0, 1, . . .). 
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and  since  limn^ooP^  =  P,  we  have  Pn  <  P, 
bounded,  therefore  it  has  accumulation  points, 
inequality  holds: 


(n  =  0, 1, . . .).  That  is,  the  sequence  Pn  is 
For  any  accumulation  point  P*  the  following 


P*  <  P.  (19) 

Let  the  sequence  Pn  defined  by  the  iteration  Pn  =  m(Pn_1)  start  from  P°  =  0.  Lemma 
1(c)  implies  m(P°)  <  m(P°),  that  is  P1  <  P1.  If  Pn  <  Pn,  then  m(Pn)  <  m(Pn)  or 
pn+i  <  p«+1.  By  induction  we  may  write 

Pn<Pn,  (n  =  0, 1, . . .). 


The  sequence  P  converges  to  P,  therefore  for  the  accumulation  points  we  have 

P*  =  P  <  P* 


(20) 


The  inequalities  (19)  and  (20)  imply  that  P*  =  P. 


□ 


V.  Joint  Power  Control,  Base  Station  Assignment,  and  Beamforming 

So  far  we  have  considered  the  power  control  problem  for  a  number  of  transmitter-receiver 
pairs  with  fixed  assignment,  which  can  be  used  in  uplink  or  downlink  in  mobile  commu¬ 
nication  systems.  In  the  uplink  power  control  problem  without  beamforming,  the  power 
allocation  and  base  station  assignment  can  be  integrated  to  attain  higher  capacity,  while 
achieving  smaller  power  allocated  to  each  mobile  as  it  has  been  demonstrated  in  previous 
studies,  [10],  [11]. 

In  the  joint  power  control  and  base  station  assignment,  a  number  of  base  stations  are 
potential  receivers  of  a  mobile  transmitter.  In  the  optimal  power  control  and  base  station 
assignment,  the  objective  is  to  determine  the  assignment  of  users  to  base  stations  which 
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minimizes  the  allocated  mobile  powers,  and  the  corresponding  transmission  powers.  Iterative 
algorithms  that  compute  the  optimal  joint  base  station  and  power  assignment  were  proposed 
in  [10],  [11]. 

In  an  uplink  scenario  where  mobiles  access  base  stations  equipped  with  array  antennas, 
the  problem  of  joint  power  control  and  beamforming  as  well  as  base  station  assignment 
naturally  arises.  A  modification  of  the  iterative  algorithm  we  proposed  for  the  joint  power 
control  and  beamforming  problem  can  perform  base  assignment  as  well.  The  iterative  step 
of  the  algorithm  is  as  follows: 

Algorithm  B 

1.  Each  base  station  in  the  allowable  set  of  a  mobile  i  minimizes  the  total  interference 
subject  to  maintaining  unity  gain  towards  the  direction  of  this  mobile: 


w, 


,n+l 


arg  mm 

w 


jjtm 


Gai(wi,a 


)Pj  +  Nm  wf w<},  (t  =  1,2, 


■jm 


m  e  Bi, 


subject  to  w/JaiTO  =  1, 


where  w”^1  is  the  optimal  beamforming  weight  vector  at  the  mth  base  station  for  the 
zth  mobile,  and  Bi  is  the  set  of  allowable  base  stations  for  the  ith  mobile. 

2.  Each  mobile  finds  the  optimal  base  station  such  that  the  allocated  power  for  the  next 
iteration  is  minimized: 


l  •  r  GkmGa  (wim  ,  a km)  n  'YoNmWim  wim  i  /„•  r>  n,r\ 

=  arg  rnm  {7o  2^ - ^ - Pk  + - ^ - />  (z  =  1,2,...,  M), 


Gin 


G, 


where  bi  is  the  optimal  assignment  for  mobile  i. 

3.  The  mobile  i  updates  the  power  based  on  the  optimum  beamforming  and  base  station 
assignment: 


pn+l  ..  VGfc6iGa,(wy1,a*6i)  „  .  _ 

Pi  =  7o  2^ - n - Pk  H - n - >  v*  —  A  A 

k^i  ^ibi 
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The  above  steps  are  combined  in  one  iteration  as 


pn+l 


mm 

w 


r  ^ km.G aiiyf imi  3-fcm)  nri  .  7o Nr 

\7o  2^ - - 1  k  + 


»w&w* 


k^i 


Gin 


G, 


(i  =  1,2,...,  M), 


subject  to  Gai( wim,aim)  =  1. 

The  convergence  and  optimality  of  this  algorithm  can  also  be  proved  following  the  same 
steps  as  in  Section  .  The  only  difference  is  that  the  minimization  is  done  over  all  beamforming 
vectors  and  base  station  assignments  instead  of  over  the  beamforming  vectors  only.  The  proof 
is  omitted  here  for  brevity. 


VI.  Simulation  Results  and  Comparisons 

We  evaluate  the  performance  of  our  algorithm  by  simulating  the  same  system  as  in  [11]. 
The  quality  constraint  is  considered  to  be  0.0304,  which  is  equivalent  to  CIR  of  -14  dB. 
This  threshold  results  in  acceptable  bit  error  rate  only  in  CDMA  systems  where  there  is  a 
processing  gain  of  the  order  of  128  or  more.  However,  the  same  methodology  can  be  applied 
to  any  wireless  network  such  as  TDMA  and  FDMA.  In  the  latter  cases  the  interference 
rejection  capability  of  antenna  arrays  can  be  utilized  to  decrease  the  reuse  distance,  or 
support  more  than  one  user  with  the  same  time  slot  or  frequency  in  each  cell.  Both  of  these 
effects  will  increase  the  capacity  significantly. 

Fig.  4  shows  a  network  with  36  base  stations  with  400  users  randomly  distributed  in  the 
area  [0.5,  6.5]  x  [0.5, 6.5]  with  uniform  distribution.  The  link  gain  is  modeled  as  Gi3  =  1  jd\^ 
where  di3  is  the  distance  between  base  i  and  mobile  j.  The  shadow  fading  is  ignored  in  these 
simulations.  Throughout  the  simulations,  we  consider  two  system  setups:  in  system  setup  I 
we  use  omnidirectional  antennas;  in  system  setup  II  we  use  antenna  array  with  4  elements. 

Fig.  4(a),  illustrates  the  use  of  system  setup  I.  Traditionally,  the  mobiles  are  assigned 
to  the  base  stations  with  the  largest  path  gains,  and  the  mobile  powers  are  obtained  by 
iterative  fixed  assignment  power  control  algorithm  as  given  by  (5).  In  Fig.  5(a),  the  dashdot 
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curve  shows  the  total  mobile  power  at  each  iteration.  This  algorithm  converges  in  about  16 
iterations.  In  Fig.  4(b),  using  the  same  system  setup,  the  base  station  assignment  is  done 
by  the  jointly  optimal  base  station  assignment  and  power  control  algorithm,  and  mobiles 
have  the  option  to  select  among  four  closest  base  stations  [11].  The  total  mobile  power  is 
depicted  in  Fig.  5(a).  The  dashed  curves  shows  that  the  total  power  is  slightly  less  than 
that  of  the  first  algorithm  considered  in  Fig.  4(a).  This  algorithm  converges  in  about  15 
iterations.  In  Fig.  4(c)  we  use  the  system  setup  II,  i.e.,  the  base  stations  are  equipped  with 
4-element  antenna  arrays.  We  apply  our  joint  power  control,  base  station  assignment,  and 
beamforming  algorithm  to  the  same  configuration  of  users  as  in  Fig.  4(a)  and  (b).  The  solid 
curve  in  Fig.  5(a)  shows  that  the  total  mobile  power  is  an  order  of  magnitude  smaller  than 
the  previous  algorithms.  Furthermore,  the  convergence  of  this  algorithm  is  much  faster;  it 
converges  in  about  5  iterations  in  our  simulation  study. 

The  capacity  of  the  system  is  defined  as  the  maximum  number  of  users  for  which  there 
exists  a  feasible  power  vector.  As  the  number  of  users  is  increased  the  maximum  eigenvalue 
of  the  gain  matrix  p(F )  approaches  unity  and  the  total  sum  of  mobile  power  is  increased. 
At  the  same  time  the  number  of  iterations  in  order  to  achieve  the  convergence  is  increased. 
Therefore,  in  our  simulations,  instead  of  setting  a  maximum  value  for  the  power  elements,  we 
set  a  maximum  of  100  for  the  number  of  iterations  required  for  convergence.  That  is,  if  the 
power  vector  does  not  converge  in  100  iterations,  we  consider  the  network  as  an  infeasible 
system. 

Using  an  antenna  array  with  four  elements  and  our  algorithm,  we  can  increase  the  capacity 
of  the  network  significantly.  In  Fig.  5(b)  the  total  mobile  power  versus  the  number  of  users 
is  depicted.  Using  omnidirectional  antenna  and  power  control  algorithm  with  fixed  base 
assignment,  we  can  tolerate  at  most  660  users.  In  the  same  configuration,  using  the  base 
station  assignment  and  power  control  algorithm  [11],  we  can  increase  the  capacity  to  800 
users.  If  we  use  antenna  array  with  four  elements,  using  our  algorithm,  the  network  can 
tolerate  2800  users.  Fig.  6  illustrates  the  base  station  assignments  for  the  above  three  cases. 
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Fig.  5(b)  shows  that  for  a  fixed  number  of  users  in  our  system,  the  total  mobile  power  is  an 
order  of  magnitude  less  than  the  other  methods. 

It  has  been  observed  in  [11]  that  integration  of  both  base  station  assignment  and  power 
control  significantly  increases  the  local  capacity,  i.e.,  handling  more  users  when  we  have  a 
hot-spot  in  the  network.  In  order  to  demonstrate  the  effectiveness  of  our  proposed  approach, 
in  Fig.  7,  400  users  are  dispersed  randomly  around  the  network.  We  then  added  users 
randomly  in  the  local  area  of  [3.5, 4.5]  x  [3.5, 4.5].  When  we  add  22  users  to  the  system 
with  setup  I,  the  traditional  fixed  base  station  assignment  reaches  its  limit.  Using  the  power 
allocation  and  base  station  assignment  [11]  and  the  same  system  setup,  when  we  add  57 
users  we  get  overload.  Using  system  setup  II  and  our  method,  we  can  add  150  users  prior 
to  overload. 

In  summary,  when  we  have  the  same  configuration  of  users,  the  use  of  4-element  antenna 
arrays  in  the  base  stations  and  our  algorithm  significantly  reduce  the  mobile  power  by  almost 
an  order  of  magnitude,  which  is  very  critical  in  terms  of  battery  life  in  mobile  sets.  It  also 
provides  faster  convergence  compared  to  the  existing  power  control  algorithms.  Furthermore, 
it  can  increase  the  capacity  of  systems  significantly. 

VII.  Conclusion 

We  have  introduced  the  consideration  of  joint  optimal  beamforming  and  power  control. 
We  provided  an  iterative  algorithm  amenable  to  distributed  implementation  which  converges 
to  the  optimal  beamforming  and  base  station  assignment  if  there  exists  at  least  one  solution 
to  the  joint  problem.  An  enhancement  of  the  algorithm  that  makes  it  appropriate  for  joint 
power  and  base  assignment  as  well  as  beamforming  was  considered  as  well. 

For  performance  evaluation  of  our  algorithm  a  notion  of  capacity  was  considered  to  be 
the  maximum  number  of  transmitters  for  which  there  exists  a  feasible  power  vector.  It  was 
shown  that  using  antenna  array  at  the  base  stations,  the  algorithm  can  improve  the  capacity 
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of  networks  significantly  in  terms  of  the  number  of  users  that  can  be  supported.  On  the 
other  hand,  it  speeds  up  the  convergence  and  saves  the  total  mobile  powers,  compared  to 
the  best  known  power  control  algorithms. 
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Fig.  1.  A  pair  of  cochannel  links  i  and  j  is  depicted. 
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Fig.  2.  Antenna  array  and  beamformer. 
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Fig.  3.  Sample  antenna  array  pattern 
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Fig.  4.  Mobile  and  base  stations  locations  for  400  users;  (a)  traditional  assignment;  (b)  optimal  base  station 
and  power  control;  (b)  optimal  base  station,  beamforming  and  power  control. 
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Fig.  5.  (a)  The  total  mobile  powers  versus  the  iteration  number;  (b)  the  total  mobile  powers  versus  the 

number  of  users. 
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Fig.  6.  Mobile  and  base  stations  locations;  (a)  traditional  assignment  with  660;  (b)  optimal  base  station 
and  power  control  with  800;  (b)  optimal  base  station,  beamforming  and  power  control  with  2800  users. 
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(a)  (b)  (c) 

Fig.  7.  Mobile  and  base  stations  locations  with  local  congested  area;  (a)  traditional  assignment  with  22 
additional  users;  (b)  optimal  base  station  and  power  control  with  57  additional  users;  (b)  optimal  base 
station,  beamforming  and  power  control  with  150  additional  users. 


